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We numerically study an anisotropic shape transformation of membranes under exter¬ 
nal forces for two-dimensional triangulated surfaces on the basis of Finsler geometry. 

The Finsler metric is defined by using a vector field, which is the tangential component 
of a three dimensional unit vector a corresponding to the tilt or some external macro¬ 
molecules on the surface of disk topology. The sigma model Hamiltonian is assumed for 
the tangential component of a with the interaction coefficient A. For large (small) A, the 
surface becomes oblong (collapsed) at relatively small bending rigidity. For the interme¬ 
diate A, the surface becomes planar. Conversely, fixing the surface with the boundary 
of area A or with the two point boundaries of distance L, we find that the variable a 
changes from random to aligned state with increasing of A or L for the intermediate re¬ 
gion of A. This implies that an internal phase transition for a is triggered not only by the 
thermal fluctuations but also by external mechanical forces. We also find that the frame 
(string) tension shows the expected scaling behavior with respect to A/N (L/N) at the 
intermediate region of A (L) where the a configuration changes between the disordered 
and ordered phases. Moreover, we find that the string tension 7 at sufficiently large A 
is considerably smaller than that at small A. This phenomenon resembles the so-called 
soft-elasticity in the liquid crystal elastomer, which is deformed by small external tensile 
forces. 

Keywords: Triangulated surfaces; Finsler geometry; Membranes; Monte Carlo simula¬ 
tions 

FACS Nos.: ll.25.-w, 64.60.-i, 68.60.-p, 87.10.-e, 87.15.ak 

1. Introduction 

A considerable number of studies have been conducted on the origins of anisotropic 
shape transformation (AST) in membranes. The AST is a typical morphological 
change in membranes, and almost all morphological changes are regarded as ASTs. 
As an example of AST, the prolate shape of fluid vesicles has its origin in the neg¬ 
ative pressure Ap inside the surface ^21^ xhe flow field is also co nside red to be the 
origin of AST for biological membranes such as red blood cells EIHSl Fluid mem¬ 
branes are expected to undergo AST when they are supported by the cytoskeleton 
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even if it is isotropic®. Also for better under stan ding of the ASTs the bilayer area 
difference elasticity (ADE) model is proposed Together with the line tension en¬ 
ergy, the ADE Hamiltonian is used for the analysis of experimentally observed ASTs 
such as the circular-to-stripe domains transition in the two-component membrane 
Moreover, the originating of the AST from the anisotropic bending ri gidity has 
been studied extensively using the Landau-Ginzburg model for membranes I ^ ^ I ^ ^ 1 IS] 

The three-dimensional structure of the lipids (tilts) or of some other molecules is 
also considered as the origin of AST The internal molecular structure influences 
the external macroscopic structure of the material. In other words, the AST is 
naturally understood as a phenomenon induced by the internal phase transition in 
the context of Einsler geometry. 

However, there is another question: whether the converse phenomena occur or 
not. The problem is whether or not the internal order/disorder transition of tilts is 
activated by AST. If this internal transition is activated, it implies that an exter¬ 
nal mechanical force can influence the internal molecular structure via the surface 
transformation, which is not only AST but also isotropic shape transformations. 
We should note that the change in the alignment of the liquid crystal molecules 
due to the external macroscopic force is considered as the origin of the so-called 
soft-elasticity in the liquid crystal elastomer (LCE), the elongation of which occurs 
with small external force ^ 1 ^ 1 ^ 

This paper aims at seeing whether or not external mechanical forces cause the 
internal order/disorder transition of cr in a Einsler geometric model, w hich is an ex¬ 
tension of Helfrich and Polyakov (HP) model for membranes I ^^ 7 1181 rppg external 
force, which is the frame tension, is supplied for fixing the surface with the bound¬ 
ary vertices 11^. The projected area Ap inside this boundary is fixed by the frame 
tension. When Ap is sufficiently large, the variable a is expected to be aligned to a 
spontaneously chosen direction. Another external force is also supplied to fix a sur¬ 
face with two vertices separated by a distance L on th e diagonal line of the surface. 
This force is identified with the string tension 7 The surface becomes oblong 
with sufficiently large L, and in this case the variable a is expected to align along 
the oblong direction. The variable a strongly aligns to this direction when the coef¬ 
ficient A of the sigma model Hamiltonian is sufficiently large, and it is expected that 
7 becomes very small in this case. This phenomenon resembles the soft-elasticity 
observed in the LCE as mentioned above. The LCE is a three-dimensional object 
and hence it is quite different from membranes, however, we expect the same mech¬ 
anism works also for the two-dimensional case. In this paper, we will confirm that 
these expected phenomena are observed in the Einsler geometric surface model. 

2. Model 

The discrete model is defined on the disk. By dividing the edges of the hexagon into 
£ pieces, we have a hexagonal lattice of size 

{N, Ne, Np) = {3f + 3^ + 1, 9f + 3£, 6f), 


( 1 ) 
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(a) 


(b) 


Fig. 1 . (a) A triangulated disk of size (A, Ap) = (91, 240,150). (b) The tangential plane at 

the vertex z, and its unit normal vector Np The variable cri(G S‘^) and its tangential component 

4 '. 

where Ne and Np are the total number of vertices, the total number of edges, 
and the total number of plaquettes (or triangles). The disk for ^ = 5 is shown in Fig. 
[TJa). The unit normal vector at the vertex i is given by 



( 2 ) 


where and are respectively the unit normal vector and the area of the 

j-th triangle linked to the vertex This uniquely defines the tangential plane 
at the vertex i (Fig. [Hb)). 



Fig. 2. (a) The triangle 123 and its three neighboring triangles, where u.i{i = 0,.., 3) are the unit 

normal vectors, and the edge vectors are given by =r 2 —ri, -^2 =i ’3 —i*i, -^3 =r 3 —r 2 . (b) The 
directions of the parameters Vij . 

The model is defined by the Hamiltonian 


S = XSo + Si + kS2, 5o =-d'’ 


cr y • cr y 



(3) 


where the variable is the tangential component of (Fig. [TJb)), and in So 
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is the sum over all nearest neighbor vertices ij linked by a bond, which is an edge 
of the triangles. We should note that ai has a value on the unit sphere G 5'^) 
and cr|' is given by 

A =ai-{ai--Ni)Ni. (4) 

Since crj' does not always belong to So in Eq. m is different from the ordinary 
sigma model Hamiltonian. and 5 ' 2 (^) in Eq* (|3|) are given by 

S,{A) = j,£l + j2£l + i3£l 

5'2(/\) = /^i(l - no • ni) + K2{1 - no • n 2 ) + /^3(1 - no • ns), (5) 


where the coefficients 7 ^ and tzi are defined by 
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The parameters Vij are given by 

Vij = 1 + [^ij] 1 ^ij — \^i ’ ijf I 5 


(6) 

(7) 


where [x] stands for Max{n G Z|n < x}, Ny = 100, and tij is the unit tangential 
vector from the vertex i to the vertex j. The inner product • tij in Eq. Q can 
also be written as cr|' • tij. We should note that Vij in Eq. © can also be written as 
'^ 7 ' = Wi • tij I + e (e = 0 . 01 ), which is almost (but not exactly) equivalent with the 
original one. This correspondence comes from the fact that Vij always appears as a 
rational form VijV^^ in the Hamiltonians. Note also that the cut off e is introduced 
to protect VijV^j^ from being divergent. The cut off e effectively makes the variable 
(Ji not perpendicular to the direction tij . 

We should note that the coefficients 7 i and Ki are defined on the edges of tri¬ 
angles. On the boundary edges the coefficients are determined only by the triangle 
in one side of the edge, while on the internal edges they are determined by the two 
neighboring triangles. Note also that Ki at the bond i becomes small (large) when 
a aligns parallel (perpendicular) to the direction of the bond i. 

The partition function is given by 

Z{X,k) = E/n drjexp[-S'((7,r)], (8) 


where the vertex number i runs from 1 to N. The center of mass of the surface is 
fixed for the free boundary surfaces. Eor the fixed boundary surfaces or the surfaces 
with fixed boundary points, i in runs only for the vertices which are not fixed. 

In the limit of A^oc, it is expected that the variable a has only the tangential 
component and hence -N^ = 0, which implies = cr|L Moreover, the magnetization 
M defined by 


M = VM • M, M = ^(Ti 


(9) 
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becomes maximal in this limit. Thus, the energy Sq effectively makes the surface 
stiffened (softened) in one direction (the other direction) at sufficiently large A, and 
hence the surface is expected to be oblong if tz is not too large. 

The discrete Hamiltonians Si and S 2 in Eq. m are obtained from the continuous 
Hamiltonians 

>51 = y ^/gdxg^^^dav ■ dbV, 

52 = i y y/g(fxg°'’’dan ■ dbn, (10) 

by replacing the metric gab with the discrete Finsler metric on a triangulated 
surface (Figs. [Ha)-[Hb)) such that^ 


F _ 
9ab 



( 11 ) 


This g^^ is obtained by deforming the diagonal elements of the Euclidean metric 
Qab = ^ab such that 1 ^ Vi 2 ^ud 1 ^ "^li^ In Si of Eq. (pd)|) . the derivatives are 
replaced by 9ir^r2—ri, 92i*^r3 —ri, where denotes the position of the vertex 
i such that >^1 = |r 2 —ri|, ^2 = |i *3 —ri| (see Figs. [2](a)-[2](b)). The derivatives in S 2 
can also be discretized in the following way: 5in->no —n 2 , 92n^no —ni, where 
rvi{i = 0,1,2,3) are the unit normal vectors shown in Fig. EJa). The integration 
measure j ^/gdPx and g^^ in the continuous Hamiltonians are respectively replaced 

by I (1/2) X;a \Jdetgh and 5“'’ (5^)“^ = where g^^ is given 

by Eq. ([TT]) . 


3. Simulation results 

Firstly, we shall describe how to update the variables in MC simulations. The vari¬ 
able cr(G S‘^) is randomly updated a a' = a Sa such that a' G 5'^, where 
the length of Sa is chosen such that the acceptance rate becomes approximately 
50% ^ 70%. The vertex position is also updated randomly r — r' + (5r so that 
the rate of acceptance is approximately 50%. One MC sweep (MCS) consists of N 
sequential updates of a and N sequential updates of X. 


3.1. Surfaces with free boundary 

In this subsection, we see the dependence of the phase structure on the parameter 
A under the fixed bending rigidity f of the model with the free boundary condition. 
As mentioned in Section El the surface is expected to be oblong (collapsed) at 
sufficiently large (small) A if is an intermediate value. Indeed, we can see that 
the surface is collapsed at A = 0.6 while it becomes oblong at A = l, and the planar 
phase appears between the collapsed and oblong phases when = 0.95 (Figs. [H^a)- 
(c)). Thus, these three different phases are expected to appear with varying A for 
f = 0.95. The phase boundary between the collapsed and planar phases is expected 
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Fig. 3. Snapshots of the surface and the surface section at (a) A = 0.6 (the collapsed phase), (b) 
A = 0.76 (the planar phase), and (c) A = 1.0 (the oblong phase). The bending rigidity is fixed to 
/c = 0.95. The red burs represent the variables CTj. The scales of the figures are the same. 


in the region 0.6 < A < 0.76, and the phase boundary between the planar and 
oblong phases is expected in the region 0.76 < A < 1. We should note that the three 
phases can be seen even when k takes the value /^ = 0.9 or k = 1. 

To see the phase boundaries more clearly, the magnetization M/N and the 
variance Cm^ which is the magnetic susceptibility, defined by 

CM = l^{{M-{M)f) (12) 

are plotted in Figs.|ll^a),(b). We find a peak in Cm at A ^ 0.86, which is the bound¬ 
ary between the planar and oblong phases. This implies that the internal transition 
of the variable a induces the shape transformation between these two phases; the 
disordered (ordered) a makes the surface planar (oblong). To the contrary, we find 
no change in the configuration of a at the phase boundary, where A 0.65, between 

the collapsed and planar phases. At A > 0.65 the magnetization M/N starts to grow 

(Fig. m^a)). The bending energy S 2 /NB and the specific heat Cs 2 defined by 

Cs, = ^((S2-(S2)f} (13) 

are plotted in Figs. |ll^c),(d). Contrary to the case of M, we see a peak in Cs 2 at 
the boundary (A 0.65) between the collapsed and planar phases. The peak is 
expected to grow larger with increasing N, and therefore the transition is expected 
to be a second order one. The configurational change in a causes this shape trans¬ 
formation transition, although the configuration of a changes continuously at this 
phase boundary. The mean square radius of gyration defined by 

i i 


(14) 
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Fig. 4. (a), (b) The magnetization M/N and its variance Cm, (c), (d) the bending energy S 2 /NB 

and the specific heat Cs 2 : cind (e), (f) the mean radius of gyration R‘^ and its variance The 

bending rigidity is fixed to /^ = 0.95. 


and the variance Cji 2 = {1/N){(^R‘^ — are plotted in Figs. |ll^e),(f). The peak 

of Cpi 2 at the phase boundary A 0.65 reflects the shape transformation transi¬ 
tion. This implies that the surface is effectively stiffened as A increases, although a 
appears to remain unchanged at A ^ 0.65 from Figs. [11(a), (b). 

Here we comment on the reason why the phase transitions are reflected only in 
certain specific quantities. The quantities S 2 and rapidly change only at A ^ 0.65, 
which is the boundary between the collapsed and planar phases, and these quantities 
change almost smoothly at A 0.86, which is the boundary between the planar and 
oblong phases, although 0^2 has a small peak at this boundary. We should note 
that the oblong phase is the one and only anisotropic phase while the other two 
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phases, the planar and collapsed phases, are isotropic. Note also that the surface 
strength, such as the surface tension and bending rigidity, along the spontaneously 
chosen direction becomes different from that along the other direction on the surface. 
From this property it follows that the anisotropic (isotropic) surface shape emerges 
as a result of the appearance of the ordered (disordered) phase of a in the context 
of Finsler geometric model. This is the reason why the magnetization M/N rapidly 
varies and hence the variance Cm has the peak only at A 0.86, which is the phase 
boundary between the oblong and planar phases. On the other hand, the fact that 
the transition at A 0.65 is reflected in S 2 and comes from the fact that the 
random configuration of a effectively softens the bending rigidity. This is another 
property of Finsler geometric model. This is possible because the surface tension 
and the bending rigidity effectively changes according to Eq. however, we don’t 
go into detail any further on this point. 

In Fig. Owe show the snapshots obtained at A = l.l for /^ = 0.95 to see the shape 
of surface in the oblong phase. We actually find that the surface is almost flat but 
it is oblong in a spontaneously chosen direction. 



(a) (b) (c) (d) 


Fig. 5. Snapshots of surfaces of size (a) A/” = 2437, (b) A/” = 3997, (c) A/” = 6487, and (c) A/” = 9577, 
obtained at A = 1.1 for n = 0.95. The variable cr, which is shown by red bur, aligns along a 
spontaneously chosen direction on the surfaces. 




0-8 ^ 



0.8 ^ 1 


Fig. 6, 


(a) The Gaussian bond potential Si/N vs. A, and (b) the triangle area A/Np vs. A. 
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Finally in this subsection, we check that 5'i/A^ = 3/2 is satisfied. This relation is 
a direct consequence of the scale invariance of the partition function for the surface 
without the boundary condition except that the center of mass of the surface is 
fixed. We find from Fig. [ 6 fa) that the expected relation Si/N = 3/2 is satisfied. 
This implies that the simulations are successfully performed. 

The mean triangle area A/Np plotted in Fig. [ 6 l^b) is slightly dependent on both 
A and N. This dependence of A/Np on A characterizes the Finsler geometric model, 
although the deviation in A/Np is very small. To the contrary, the triangle area 
is always constant in the canonical surface model, which is defined by Si = 
and S 2 = ^^j{l — ni • n^). In this canonical surface model, the relation S'i/A^ = 3/2 
automatically implies that A/Np remains constant, because the regular triangles 
always dominate the ensemble configurations. 

An important feature of the model considered in this paper - that A/Np is not 
constant in the model of this paper - comes from the fact that the surface tension 
coefficients 7 ^, (i = 1,2,3) in Eq. (j 6 j) are not constant and vary depending on the 
position and the direction on the surface. As a consequence, the bond length squares 

(i = l, 2, 3), and hence the triangle area becomes dependent on the position and 
the direction on the surface. The dependence of A/Np on N indicates that 7 ^ or Vij 
in Eq. © are strongly influenced by the size effect of the surface size N. Indeed, 
the linear plot of A/Np vs. 1/N at A ^ 0.7 shown inside Eig. [61(b) indicates that 
A/Np has a finite value in the limit of N ^oc. This size effect implies that the long 
wave length modes in the fluctuations of a are important for the correlation energy 
between in 5'o- 

3.2. Frame tension of surfaces with fixed boundary 

In this subsection, we see how the shape transformation is reflected in the internal 
degree of freedom a. To see this, all boundary vertices are fixed on a plane so that 
the projected area Ap enclosed by the boundary becomes a constant Ap{= aNp), 
where 


a = Ap/Np (15) 

is the area of the projected triangle. It doesn’t depend on the lattice size N. This 
implies that the length of the diagonal line of the lattice of size N^ such as shown 
in Fig. HJa), is given by L = 2{a /where i is the partition number of the 
hexagon edge introduced at the beginning of Section |2l This constant area a is 
varied while A and k are fixed, and the dependence of the a configuration on a is 
checked. We should note that Ap is not always identical with the real surface area, 
which varies due to the surface fluctuations. 

The boundary condition assumed in this subsection is used to calculate the 
frame tension r. As a consequence, we will see how the external force r influences 
the internal variable a. 

Eigures [3(a),(b) show snapshots of the surfaces of size = 9577 for a 0.027 
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Fig. 7. Snapshot of the surfaces oi N = 9577 for (a) a ~ 0.027 and (b) a ~ 0.086, where A = 0.9 
and Av = 0.65. The variable a is not shown on the surfaces. 


and a 0.086, which respectively correspond to the boundary hexagon diameter 
Lc^lS.S and L^^SO, obtained under A = 0.9 and /^ = 0.65. 



Fig. 8 . (a) The magnetization M/N vs. the projected triangle area a, (b) the variance Cm vs. a, 

(c) the absolute internal energy |S'o|/A/' vs. a{N/No)^-^, and (d) the frame tension (r—ro){N/No)^-^ 
vs. a{N/NQ)^-^^ , where to = —0.2 and Aro = 2437. 


The magnetization M/N and its variance Cm are plotted in Figs. [Hfa),(b), and 
we find that the peak of Cm increases with increasing N. This implies that the 
external force, which fixes the surface boundary, can influence the internal variable 
cr. Moreover, we And that M/N scales as a function of a at the transition region, 
which is shown in Fig. [Sfa). Apparently, this is a nontrivial result, because there 
is exactly no dependence of M on the projected area in the case of the ordinal 
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sigma model, which is a model of the ferro-magnetic transition. We also expect 
from Figs. [Hl^a),(b) that a is ordered (disordered) in the limit of a ^ oo (a ^ 0). 
The rotational symmetry of a is spontaneously broken if the surface is expanded to 
have sufficiently large area. This is in sharp contrast to the original sigma model, 
in which the symmetry is spontaneously broken only when A ^ oo. The internal 
energy Sq/N is also seen to scale as a function of a{N/No)^'^ as shown in Fig.[ 8 l^c). 
The frame tension r is defined as in MM 

T = {2{Si) - 3N) / {2Ap). (16) 

This relation comes from the scale invariance of the partition function, which is 
given by dalog Z{aT)\a=i =0. Indeed, using the fact that S 2 is scale independent 
and 5'i((ar) = (a^5'i(r), and that da \og Z{Ap; (ar)|cK=i =da \og Z{a~^ Ap] r)|cK=i can 
be written as —2ApdAp log Z{Ap; r), and the assumption that the partition function 
for the macroscopic membranes is given by Zmac = exp {—F) with the free energy 
F = tAp^ we have r in Eq. (pffi|) xhe results in Fig.[Hfd), shown in a linear scale, 
indicate that {r — tq){N with tq = — 0.2 scales as a function of a{N 
at the transition region. 

3.3. String tension of surfaces with two-point boundaries 


(a) (b) 



Fig. 9. Snapshot of the surfaces of A/” = 6487 for (a) L = 40 and (b) L = 80, where the parameters 
are fixed to A = 0.84, /^ = 0.75. The variable a is not shown on the surfaces. 


The calculation for the string tension measurements needs the boundary condi¬ 
tion for the surface, which is of a hexagon, such that the two boundary vertices are 
fixed with a sufficiently large separation distance. The two vertices are on a diago¬ 
nal line of length L. This implies that the string tension 7 is applied to the surface 
as an external force to make the surface oblong so that the distance between the 
vertices is fixed to L. Snapshots in Figs.[9l^a), (b) obtained with I/ = 40 and L = S0 
are relatively oblong, although the central parts are not so thin. 

We should note that almost all triangles become oblong on the surface but small 
part of triangles or vertices, which are located far from the diagonal line between the 
two terminal points, are influenced only slightly by the extension. For this reason. 
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the scaling behavior of the string tension a with respect to N is expected to deviate 
from the ordinal one 123^ 



Fig. 10. (a) The magnetization M/N vs. L/-\/N, (b) the variance Cm vs. L/-\/N, (c) the internal 

energy |S'o|/A/’ vs. and (d) a log-log plot of the string tension 7 vs. where 

7 ~ is expected. 

The magnetization M/N and the variance Cm vs. L/\/N are plotted in Figs. 
[TOlf a). (b). The reason why L/\/N is used is because L/\/N plays a role of the pro¬ 
jected chain length IlH the oblong surface is identified with a linear chain of length 
L = y/NiL- This L/y/N) is expected to be independent of N. It corresponds 
to the projected triangle area a introduced in Section [3^ for the calculation of the 
frame tension. The peak of Cm increases with increasing N (Fig. flOl b)). and this 
implies that the external force causes an internal transition of a as in the case for 
the frame tension in Section IS21 Thus, the symmetry for a is also spontaneously 
broken if the surface is expanded to be oblong by external forces. 

Now the string tension 7 is again defined as in |20l [2T] 

7 = (2(5i)-3V)/(3L). (17) 

This formula for 7 is obtained by the scale invariance of Z for the surfaces with 
fixed two point boundaries separated by L. 

The internal energy |5'o|/A^ and the string tension 7 are expected to scale as a 
function of L/N^. Indeed, we find that |5'o|/A^ and 7 scale as a function of 
and respectively as shown in Figs. flOlc). (d). We find that the string 

tension, plotted in a log-log scale in Fig. flOld). scales as 7 ~ although 
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the range of the variable is not sufficiently large. This scaling behavior 

deviates from the ordinary one such as 7 ^ {L/NY (see [20]). We suppose that the 
reason of this deviation is that a small number of vertices does not always contribute 
to the string tension r as mentioned above. 


3.4. String tension at high and low A 



Fig. 11. (a) The string tension 7 vs. L/\fN obtained with A = 0 and A = 5, and (b) log-log plots 

for a scaling behavior of 7 , where A = 0 and A = 5. The bending rigidity is fixed to /c = 0.25 in both 
(a) and (b). 


From the results presented in Figs. |4]^a)-(f), we have found that the surface 
becomes oblong at sufficiently large A for the small k region. In this case, the 
variables a align in the oblong direction. This indicates that the string tension 7 , 
which is identified with the external force applied to the surface to maintain the 
surface oblong, is influenced by the presence of the variable a at sufficiently large 
A. 

In this subsection, we will show that there is a nontrivial difference between 
7 ’s at large and small A. Indeed, 7 vs. L/\/N plotted in Fig. fTTIa) has a large 
difference in their values, where A = 0 and A = 5. The bending rigidity is fixed to 
a relatively small value n = ^.2b. The difference of these two 7 ’s remains constant 
and independent of L/\/N at least for the region 0.5 <I//a/]V <1.5. Therefore, we 
understand that the difference is quite large compared to 7 itself at least for the 
small L/\YN region (see Fig. fTTIa)). To see this difference more clearly, we plot 
7 vs. in a log-log scale in Fig. fTTIb). We see that the difference becomes 

zero in the sufficiently large ^ oc. To the contrary, we also see that the 

difference becomes larger for the smaller . This implies that the surface 

becomes oblong with small tensile force for the large A region. This phenomenon 
resembles the soft-elasticity seen in the LCE g^j^d therefore we consider that 
the Finsler geometric model for membranes can be extended to a model for three 
dimensional objects such as the LCE. 
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4. Summary and conclusion 

We have studied a Finsler geometric model to find how the external mechanical 
forces influence the internal variable a, which represents the tilt and is used to 
define the Finsler metric. The sigma model interaction energy for the tangential 
component of a is included into the Hamiltonian with the interaction coefficient A. 

Firstly, we find that the model undergoes an anisotropic shape transformation 
when A changes in the range 0.5< A< 1.5 under a fixed bending rigidity k which is 
relatively small. The surface becomes oblong in the limit of A^oo; but it becomes 
planar and then collapsed as A reduces to A^O. Secondly, we find that the frame 
tension r, as an external force, causes the internal transition of a such that a 
becomes ordered (disordered) when r is sufficiently large (small), which corresponds 
to large (small) projected triangle area a. Thirdly, it is also found that the string 
tension, as an external force, causes the internal transition of a such that this vector 
becomes ordered (disordered) when the distance L increases (decreases) between 
the boundaries. Thus the numerical results presented in this paper confirm that 
the internal order/disorder transition causes the shape transformations and vice 
versa in the context of Finsler geometric surface model for membranes. Moreover, 
the string tension 7 at the large A region is considerably smaller than 7 at A—>0. 
This phenomenon is expected to share the same origin with the soft-elasticity in 
the liquid crystal elastomer (LCE), and therefore this indicates the possibility that 
the Finsler geometric model for membranes is extended to a model for 3D objects 
such as the LCE. 
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